A 173  477 


NPS-53-86-01 2 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


DTIC 

cfc‘£LECT? 

i*A  OCT  2  7  1988 


OCT  2  7  1986 


SPECIAL  METHODS  FOR  PROBLEMS  WHOSE 
OSCILLATORY  SOLUTION  IS  DAMPED 


Beny  Neta 
September  1986 
Technical  Report  For  Period 
April  1986  -  September  1986 


Approved  for  public  release;  distribution  unlimited 

Prepared  for:  Naval  Postgraduate  School 
Monterey,  CA  93943 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CALIFORNIA  93943 


R.  C.  AUSTIN 

Rear  Admiral,  U.  S.  Navy 
Superintendent 


D.  A.  SCHRADY 
Provost 


-Reproduction  of  all  or  part  of  this  report  is  authorized. 
This  report  was  prepared  by: 


r\Jt 


BENY  NETA 

Associate  Professor  of 
Mathematics 


Reviewed  by: 


Released  by: 


M.  FREDRICKSEN,  Chairman 
Department  of  Mathematics 


RSHALL 
Dean  of  Information 
and  Policy  Sciences 


<!< 

•!* 

I 


$ 

ft 

C», 

V 


lY 

lY. 


m 

$ 

w 

w 

w 

•s* 

»V* 

w 

S' 

h 


m 

mmmm. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RACE  (Whit  Data  Enl.r.tfJ 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

NPS-53-86-01 2 


4.  TITLE  f«id  Submit) 

Special  Methods  for  Problems  Whose 
Oscillatory  Solution  is  Damped 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


S.  TYRE  OF  REPORT  A  RERIOO  COVERED 

Technical  Report 
04/86  -  09/86 


4.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AuTMORfa; 


Beny  Neta 


•  ■  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA  93943 


II.  CONTROLLING  office  name  AND  AOORESS 


ontract  or  grant  numbers 


i 


12.  REPORT  DATE 

September  1986 


IS.  NUMBER  OF  PAGES 

19 


MONITORING  AGENCY  NAME  *  ADDRESS (II  dttloront  from  Controlling  Olllco)  IS.  SECURITY  CLASS,  (ol  thlo  roport) 

.  Unclassified 


o.  DECLASSIFICATION/ DOWNGRADING 
SCHEOULE 


IS.  OlSTRI EUTlON  STATEMENT  (ot  thlo  Roport) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (ol  tho  obotroct  entered  In  Block  20,  II  dllloront  from  Roport) 


19.  KEY  WORDS  (Conllmi*  on  rtr«r««  oldo  II  nocoooory  ond  Idontlfy  by  block  number) 

ordinary  differential  equations,  exponentially-fitted  methods,  shallow  water 
equations,  damped  oscillations 


20.  ABSTRACT  fConllnu#  on  rereroe  oldo  If  nocoooory  ond  Identify  oy  block  number) 

■  This  paper  introduces  methods  tailored  especially  for  problems  whose 

solution  behave  like  <0$  ,  where  i  is  complex.  The  shallow  water 
equations  with  topography  admit  such  solution. 

This  paper  complements  the  results  of  Pratt  and  others  on  exponential - 
fitted  methods  and  those  of  Gautschi,  Neta,  van  der  Houwen  and  others  on 
trigonometrically-fitted  methods. 


DO  ,  j  an*7I  1473  EOITION  OF  I  NOV  •«  IS  OBSOLETE 


S  N  0102-  LF-014-6401 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  of  THIS  PAOE  (Whu  Omtm  tntimd) 


ABSTRACT 


This  paper  introduces  methods  tailored  especially  for  problems 
whose  solution  behaves  like  e  x,  where  \  is  complex.  The  shallow 
water  equations  with  topography  admit  such  solution. 

This  paper  complements  the  results  of  Pratt  and  others  on 
exponential-fitted  methods  and  those  of  Gautschi,  Neta,  van  der 
Houwen  and  others  on  trigonometrically-fitted  methods. 

1 .  Introduction 

In  this  paper  we  consider  linear  multistep  methods 


a.y  ,  =  h 

;  =  0  *•  n+1-x-  ,>o 


n+i-. 


.  )  ,  k  *•  1,  n  >  k-1 


for  integrating  the  initial  value  problem 


y  '  ( x )  =  f  (x ,  y  (x )  )  ,  y(xn)  =  yn  . 


?his  linear  multistep  method  is  characterized  by  the  polynomials 


^  V  _  n 

( , )  =  >  a.c  ,  c 

>  =  0 


(')  = 


2  b. 

:.  =  0 


The  main  assumption  of  this  paper  is  that  it  is  a  priori  known 
that  the  solution  is  approximately  of  the  form 


y  (x)  -  cQ  + 


i>  .  t 
.  1 


0U»t'TY 
IN»R  TEO 


4= 
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'  i'/j; id,..  .  _  „2(fCS 
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where  =  wj  +  j •  and  the  frequencies  w^  are  in  a  given  interval 
[w£,wu]. 

The  special  case  where  A^  =  jwg  with  w^  given  was  considered 
first  by  Gautschi  [8].  His  approach  was  the  following.  Let: 


(2 )  =  c  ( ez )  -  zc  (ez ) 


(5) 


then  the  local  truncation  error  of  (1)  is  given  by  Lambert  [11] 


Ln+k 


:(h  y(tn) 


(6) 


Insertina  (4)  in  (6)  vields 


Tn+k  '  :  (0)c0 


m  i  A  .  t 

~  c  .  t •  ( ih  A  . )  e  ^ 

1  =  1  ^ 


'  Xj  =  jv0 


;7) 


The  coefficients  b^  are  chosen  in  such  a  way  that 

: (ihjw0)  =  0  ,  j  =  0,1,. ..,q  ,  (8) 

for  the  largest  value  of  q  possible.  q  is  then  called  the  trigo¬ 
nometric  order  of  the  method.  Gautschi  has  chosen  a^^  such  that  the 
methods  are  of  Adams  and  Stormer  type.  However,  these  methods  are 
sensitive. to  changes  in  the  frequency  Wq .  Neta  and  Ford  [13] 
developed  Nystrom  and  generalized  Milne-Simpson  type  methods.  These 
methods  showed  less  sensitivity  to  perturbation  in  Wq  but  require 
the  eigenvalues  of  the  Jacobian  to  be  purely  imaainary.  Neta  [14] 
has  developed  families  of  backward  differentiation  methods  that 
overcome  the  above-mentioned  restriction.  Salzer  [17]  has  developed 


predictor-corrector  methods  based  on  trigonometric  polynomials. 
See  also  Steifel  and  Bettis  [18]  and  Bettis  [3] .  Van  der  Houwen 
and  Sommeijer  [10]  have  developed  an  alternative  approach.  The 
conditions  (8)  were  replaced  by 


*  (0)  =  0  , 

(9) 


<t  (ihX  (j)  )  =  0  ,  j  =  1,2,.  .  .  ,q 


where  the  X  ^ ^  are  appropriately  chosen  points  in  the  interval 
[w^wu]  . 

An  advantage  of  this  so-called  minimax  approach  over  the  fitting 
approach  is  the  increased  accuracy  in  cases  where  no  accurate  esti¬ 
mate  of  w q  is  available  or  when  the  frequency  is  varying  in  time. 

Tne  ether  special  case  considered  in  the  literature  is  where 
Xj  =  i: y  Probably  the  first  article  on  the  subject  is  due  to  Brock 
and  Murray  [5].  They  discuss  the  use  of  exponential  sums  in  the 
integration  of  a  system,  of  first  order  ordinary  differential  equations 
Dennis  [7]  also  suggested  special  methods  for  problems  whose  solution 
is  exponential.  He  suggested  a  transformation  of  variables.  More 
recently,  Carroll  [6]  has  developed  exponentially  fitted  one-step 
methods  for  the  scalar  Riccati  equation.  For  the  general  first 
order  system  of  equations,  Pratt  [16]  suggests  methods  based  on  the 
three  parameter  exponential  function 


I  (x)  =  A  +  Be 


zx 


(10) 


The  parameters  A,  B  are  given  in  terms  of  values  of  y  and  f. 

Several  possibilities  for  z  are  given  based  on  results  of  Brandon 
[2]  and  Babcock  et  al .  [1J . 

Lyche  [12]  analyzes  multistep  methods  which  exactly  integrate 
w  x 

the  set  {x  e  },  where  wn  is  real  or  imaginary. 

In  this  article  we  developed  various  methods  fitting  exponentials 
and  methods  obtained  via  the  minimax  approach. 

2 .  Construction  of  Methods 
2 . 1  Fitting  Methods 

In  this  subsection  we  discuss  various  fitting  methods.  To  this 
end,  we  separate  MihjX)  =  0,  j  =  l,2,...,q,  into  real  and  imaginary 
parts.  This  yields  the  following  equations  relating  the  coefficients 
a  -  ,  b  .  , 


a  .  e- 


(k- 


=  0 


cos  j v  (k- ^ ) 


l  a.e^"^  '  *  sm  jv(k-I) 

x  =  0  '■ 


J  b,e^'  (K“V }  [  j u  cos  jv  (k-1) 
:=0  " 


-  jv  sin  j v  (k-i ) ] 


7  b„e-^  ^  1  ^  [  ju  sin  jv  (k-i  ) 

;.=o  * 

+  jv  cos  jv  (k-i) ] 


0  , 

(11) 


0  , 


where  =  w  +  i^,  u  =  -h-;  ,  v  =  hw,  j  =  1,2,...,  a. 

For  explicit  methods,  bQ  =  0.  For  Adams  type  methods  aQ  =  1,  a^.  =  -1, 
a^  =  0  for  i  =  2,...,k.  For  h’ystrom  or  generalized  Milne-3impson 
methods  aQ  =  1,  a2  =  -1  and  other  a^  = 


0. 


k  =  1  Implicit 


Adams 


4< 

m 


viV 

#4 

Vt.V 


w 

v,vi 

,fS! 


M 

m 


v;>1 


>4 


ve  +  y  smv  -  v  cosv 
2  2 

(y  +vi")  sin  v 


ve  -  u  smv  -  v  cosv 


(y  +v  )  sin  v 


1  ~ '  COS  V 

For  -  =  0,  the  coefficients  become  bn  =  - > -  =  b,  ,  which  aaree 

0  sm  v  1' 

with  Gautschi  [8]  if  the  coefficients  are  expanded  in  Taylor  series 
with  respect  to  v. 

k  =  2  Explicit 
Adams 


= 


(t  sin2v  -  v  cos2v)e"  +  v  cosv  -  t  sinv 

,2  2, 

(t  +.  )  sm  v 


( v  cos  v  -  t  sin  v  )  e 
,  2^  2. 

(y  +-.  )  sm  v 


Kv strom 


(y  sin2v  -  v  cos2v)e" 
,  2  2, 

(y  +v  )  sm  v 


,  ,  2  y  -  y 

,  (v  cos v  -  u  sin v) e  +  ve 

°2  ~  - 2 — T - 

(y  +v  )  sin  v 


s  in 

For  t  =  0,  the  coefficients  become  b^  =  2  — - — b^  =  0  which 


aqree  with  Neta  [14] 
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k  =  2  Implicit 

In  this  case,  one  obtains  a  one-parameter  family  of  (Adams, 
generalized  Milne-Simpson)  methods  of  trigonometric  order  1.  The 
free  parameter  can  be  used  to  increase  the  algebraic  order  of  the 
method  as  in  [1 3] . 

Backward  Differentiation 

e^u  cos  2v  +  a^eu  cos  v  +  -  b^e^^tu  cos2v  -  v  sin2v)  =  0  , 


2 

e  "  sin  2v  +  a^e"  sin  v 


-  b^e  “(u  sin2v  +  v  cos2\)  =  0  ,  (15) 


=  0  . 


This  system  can  be  solved  by  MACSYMA  (Project  MAC'S  SYmbolic  MAnipu- 
lation  system  written  in  LISP  and  used  for  performing  symbolic  as 
well  as  numerical  mathematical  manipulation  [4])  or  by  REDUCE  [9]. 
The  solution  is 


al  = 


a2  = 


_ ve  -  u  sin  2-;  -v  cos2v _ 

-e"  (v  cosv  +  p  sinv)  +  p  sin  2  .  +  v  cos2v 


•1  “  aL  , 


.  -e  sin  \i  +  e“  sin2v  -  sin  v 

b  =  — T - ; - ; - 

-e  p(v  cosv  +  p  sinv)  +  eM  (p  sin2v  +  v  cos2v) 


For  1  =  0,  the  coefficients  agree  with  those  given  in  [14] 


Again  here,  one  obtains  a  one-parameter  family  of  methods  of 
trigonometric  order  1.  In  order  to  get  methods  of  trigonometric 
order  2,  one  has  to  construct  a  3  step  implicit  method  of  Adams  or 
generalized  Milne-Simpson  type.  In  order  to  increase  the  trigo¬ 
nometric  order  without  going  to  a  higher  step  number,  one  can  con¬ 
struct  linear  multistep  methods  for  which  the  coefficients  a „  are 
also  functions  of  ; ,  w.  Some  examples  are  given  in  the  next 
subsection . 

2 . 2  Generalized  Fitting  Methods 

In  this  section,  we  construct  some  linear  mulstistep  methods 
of  the  form 

;-0  a;  (/')yn+l-;,  =  h  £|0b£  {A)  fn-fl-£  '  k  -  1'  n  -  k_1  •  (17) 

Since  a.  are  functions  of  X  one  has  more  free  parameters  for  his 
disposal  which  can  be  used  to  obtain  higher  trigonometric  order 
methods  with  relatively  lower  step  number. 

k  =  2  Implicit 

In  this  case,  one  has  to  solve  the  following  linear  system  of 
five  equations  for  the  parameters  a^,  a2 ,  bg ,  b^ ,  b2  to  obtain  a 
method  of  trigonometric  order  2. 


0 


1  +  a , 


cos2v  +a^eH  cos  v  +a2  -bge^(u  cos2v  -v  sin2v) 

-b^eu(u  cosv  -v  sinv)  -  ub2  =  0  , 

2 ,  ,  2  , 

e  "  sin2v  +a^G4  sinv  -b^e^  (y  sin2v  +v  cos2v) 

-b^e^(u  sinv  +v  cosv)  -vb2  =  0  , 

4 '  ■  2  u  4 1 1 

e  “  cost.  +a^e  cos2v+a2~bQe  "  (2u  cos  4v  -  2v  sin4v) 

2’- 

-b^e  “  (2y  cos2v  -2v  sin2v)-2yb2  =0  , 

4-  2u  4> 

e  “  sin4-.  -|-a.|e  sin2v  ~^oe  (2u  sin4v  +2v  cos4v) 

-b^e  "{2,.  sir.2v  +2v  cos2v)-2vb2  =  0. 


The  system  was  solved  by  REDUCE  [9].  The  expressions  for  the 
coefficients  are  complicated  but  REDUCE  produces  an  output  in  the 
form  of  Fortran  statements  that  can  be  incorporated  in  a  computer 
crccram.  for  numerical  experiments  with  such  a  method. 


2  .  3  y.ir.imax  Methods 

In  this  section  we  discuss  minimax  methods,  i.e.,  methods 
obtained  by  satisfying  conditions  (9).  These  conditions  can  be 
written  in  terms  of  a„,  b„  as  follow’s: 


(k-I  )  .. 


cos  (k-1 ) v  ^ 


-  b.e(k-^  J^-v(^sin(k-,)v(^+Jj)cos(k-;)v 

,=0  ' 

j  =  1,2,. . . ,q  , 


_  -  i  ( 3  >  a 


.  y  b  .U-:)s<5,{vtJ)co#u.t)v<j)+Bc)»,lo(lt.{)v(j)> 

£  =  0 


3  r 


where  t*’1  =  h-;(3> 


v  <3  >  -  hw*!* 


ihAVJ  are  the  zeros  of  the  function  c(ih')  such  that  it  has 


a  small  maximum  norm  in  the  rectangle  £  w  _<  wu,  v/  _<  v  _< 


To  obtain  the  best  approximation  in  this  case  is  certainly 


not  easy;  but  we  will  assume  that  one  can  write  <j>  (ih\  ^  )  = 


as  a  product  of  2  one  variable  functions.  Thus 


and  w(3}  can  be  taken  as  Chebyshev's  points  on  the  corres¬ 


ponding  interval,  i.e. 


(j)  =  k(-  +•-..)+  i($  *Cp)  cos 


\  -  .  ’  -  .  .  /  '  o'r'.  V/7, 

2  l  u  2  u  't 


(j)  1.  ,  1,  .  2j-l 

W  =  2(W1  +Wu]  +  2{wu  -w£)cos  -fe" 


1  —  1,2,... ,g  . 


For  this  choice  of  points,  one  can  evaluate  the  coefficients 


a  n  ,  b ,,  by  solving 


C  ( 0 )  =  0  , 


,(r(^,w(j))  =  0  ,  j  =  1,2,... ,q 


We  call  such  methods  product  minimax  (PM  ) . 


The  number  of  free  parameters  for  implicit  methods  2k+l  and 
the  number  of  equations  is  2q+l,  thus,  the  trigonometric  order 
q  is  equal  to  the  step  number  k. 

k  =  1  Implicit 
In  this  case 


*(1)  =  +v 

w(1)  =  +wu)  (23) 


a 


0 


1  = 


-a 


1 


and  the  system  of  equations  can  be  solved  for  b^,  b^.  This 
yields  the  coefficients  given  by  (12)  where 


u  =  -h; 


(1) 


v  -  hw 


:d 


(24) 


Thus,  the  product  minimax  method  would  suggest  using  the  center 

of  the  rectangle  Iv^,^u)  x  [w^,  ,wu)  as  Aq.  To  obtain  a  product 

minimax  method  of  trigonometric  order  2,  one  has  to  solve  a 

system  of  5  equations  similar  to  equation  (18)  with  the  unknowns 

bQ,  b^,  b2 ,  a^,  a2 .  The  difference  is  that  in  the  last  2  equa- 

(2)  (2) 

tions  one  should  replace  2p  by  p  and  2v  bv  v  .  In  the 

second  and  third  equations  of  (18),  the  p,  v,  should  be  replaced 
by  p(1) ,  v{1)  respectively.  The  resulting  system  can  be  solved 
by  MACSYMA  [4]  or  REDUCE  [9], 

In  the  next  section  we  implement  two  methods  of  trigonometric 
order  1  and  2  and  see  how  the  product  minimax  methods  compare 
with  fitting  methods. 


11 


3 .  Numerical  Example 

Both  systems  (18)  and  (19) -(21)  were  solved  by  REDUCE  which 
produced  a  FORTRAN  subroutine  for  the  evaluation  of  the  coeffi¬ 
cients.  This  subroutine  is  called  only  once  during  the  integration. 

The  methods  were  compared  for  the  solution  of  the  initial 
value  problem 


z  -  j  i  ( 1  +  i )  z  =  0  ,  0  £  t  <_  4  , 

z(0)  =  1 


(25) 


whose  exact  solution 


iy(l+i) t 

z (t)  =  e 


thus 


(26) 


=  /  =  ~  J  >  w  =  j  .  (27) 

In  order  to  avoid  complex  arithmetic,  we  rewrite  the  differential 
equation  as  a  system  of  equations  for  the  real  and  imaginary 
part  of  z  =  u  +  iv. 


u  +  ( u  +v)  =  0  , 
v  -  ^  ( u  -v)  =  0  , 


v (0 )  =  0  . 


0  <  t  <  4 


(28) 
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H 


The  system  is  solved  by  fitting  methods  of  trigonometric  order 
1  and  2  with  h  =  .  01  and  various  values  of  ip  and  w.  In  Table 
1  we  list  the  Euclidean  norm  of  the  error  at  t  =  4.  It  is  • 

clear  that  the  method  is  not  sensitive  to  perturbations  in  the 
values  of  tp  and  w. 


Ail' 

L  w 

error 

first  order 

second  order 

1 

0 

0 

. 3678 (-12) 

.4433 (-12) 

i  0 

.1 

.4482  (-7) 

.3508 (-11) 

!  -1 

0 

.4346 (-7) 

.2544  (-11) 

!  .i 

i 

.  1 

.6342 (-7) 

.  4421  (-11 ) 

0 

.2 

.9241 (-7) 

. 8126  (-11 ) 

.2 

0 

.8706  (-7) 

.5095 (-11) 

Table  1 


Using  the  product  minimax  methods  of  trigonometric  order 
1  and  2  with  h  =  .01  and  various  squares  centered  at  ^  =  -  tt/2, 
w  =  tt/2,  the  error  is  much  laraer  but  again  is  insensitive 
to  small  perturbations  in  the  length  of  the  sides  of  the 
squares.  In  Table  2  we  list  the  Euclidean  norm  of  the  error 
at  t  =  4. 
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aMMMMmililiil . . 


length  of  side 


first  order 


second  order 


.4 
.  8 
1.2 
1.6 
2 


3678  (-12) 
3678  (-12) 
3678  (-12) 
3678 (-12) 
3678  (-12) 


.1469 (-6) 
.2348  (-5) 
.1186  (-4) 
.  3728  (-4) 
.9001  (-4) 


Table  2 


Note  that  the  perturbations  in  the  product  minimax  methods  are 

larger  than  those  allowed  in  the  fitting  methods.  It  is 

possible  that  the  larger  errors  in  the  product  minimax  methods 

are  due  to  the  assumption  that  ?  can  be  written  as  a  product 

of  2  one  variable  functions.  Also  note  that  for  the  first 

2 

order  method,  one  always  gets  a  good  result  since  PM  always 
uses  the  center  of  the  square. 
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APPENDIX 


Here  we  show  that  the  shallow  water  equations  with 

A  x 

topography  have  a  solution  of  the  form  e  ,  where  A  is  complex. 
This  system  of  equations  consists  of  three  equations  with 
three  forecast  variables  u,  v  and  <f .  The  equations  are: 


Su  J  Su  ,  Su  ,  ,  Sf  n 

tt—  +  U  ■= —  +  V  -r—  -  fv  +  -  0  , 

St  dX  Sy  ' 


dX 


|i  +  u|i.viv  +  fu  +  ii=o, 

ot  Sx  ay  ay 


i  +  +  klvl*':e)l  -  0  • 


( A .  1 ) 


(A. 2) 


(A. 3) 


where  <t  =  gh  is  the  geopotential  height  (h  =  height  of  free 
surface) ,  is  the  bottom  topography  (assumed  to  be  independent 
of  time),  u,  v  are  the  components  of  the  wind  velocity  in  the 
x,  y  direction,  respectively,  and  f  is  the  Coriolis  parameter. 
Linearizing  the  equations  by  letting 


u  =  U  +  u '  ,  v  =  V  +  v '  ,  ;  =  4>+  i  '  , 

where  U,  V  are  the  constant  mean  flow  and  i>  is  independent  of 
time.  Assuming  that  U,  V  are  related  to  <P  via  the  geostrophic 
relations 


u  -  .  k  21 

U  f  dy  ' 


v  =  l  £* 

v  f  Sx 


(A. 4) 


one  obtains  the  linear  svstem  (after  dropping  the  primes) : 
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where  y  =  $  -  <f 

D 

If  the  flow  is  assumed  to  be  along  the  topography  as  in 
[15],  then  the  right  hand  side  of  (A. 7)  is  zero.  In  such  a 
case,  one  can  write  the  solution  in  the  form 


u 


v 


u  e1  (Cx  +,~y  "Gt) 
o 

i  ( £x  +  nv  -  ct ) 
v  e  - 
o 

c  ei  +r,y  -ct) 
To 


where 

S  =  w  ~  ip  , 
n  =  v  -  ie  . 

In  order  for  (A. 8)  to  be  a  solution  for  (A.5)-(A.7),  one  must 
have 


A  =  -o  +  £u  +  nv 


satisfying 


(A. 5) 

(A. 6) 

(A. 7) 


(A. 8) 

(A.  9) 

(A. 10) 
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i  A  [  ( i  A )  2  -  in(inY  +  |^-1  -  f  t  -  i  A  f  -  (inY  +  |y)  ] 

+  (i£Y  +  |^}  ( -in f  +  Z,\)  =  0  . 


The  real  and  imaginary  parts  of 


>  =  -a  +  yU  +  vV  , 


\i  =  -cU  -  8 V  , 


(A. 11) 


(A. 12) 


satisfy  the  following  system  of  equations  (after  dropping 
nonlinear  terms  in  A) 


>r[f*+vU2+n2>l  *  &>  ■  li-  «  • 


«■:  -  >i«2  +v<;2  ^2»i  *  “  • 


(A. 13) 


In  general,  "■  is  complex  and,  thus,  the  shallow  water  equations 
have  a  solution  in  the  class  of  problems  to  be  discussed  here. 
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